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Abstract 

Braneworld models with variable brane tension A introduce a new degree of freedom that allows 
for evolving gravitational and cosmological constants, the latter being a natural candidate for dark 
energy. We consider a thermodynamic interpretation of the varying brane tension models, by 
^ ' showing that the field equations with variable A can be interpreted as describing matter creation 

^H \ 

^D ■ in a cosmological framework. The particle creation rate is determined by the variation rate of the 

^^ , brane tension, as well as by the brane-bulk energy-matter transfer rate. We investigate the effect 

vn 

<^ ' of a variable brane tension on the cosmological evolution of the Universe, in the framework of a 

o ■ 

particular model in which the brane tension is an exponentially dependent function of the scale 
• >— 1 , factor. The resulting cosmology shows the presence of an initial inflationary expansion, followed 

X, 

H by a decelerating phase, and by a smooth transition towards a late accelerated de Sitter type 

expansion. The varying brane tension is also responsible for the generation of the matter in the 
Universe (reheating period). The physical constraints on the model parameters, resulted from the 
observational cosmological data, are also investigated. 
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I. INTRODUCTION 

The idea of embedding our Universe in a higher dimensional space has attracted a con 
siderable interest recently, due to the proposal by Randall and Sundrum that our four- 



Q&. 



dimensional (4D) spacetime is a three-brane, embedded in a 5D spacetime (the bulk) 
This proposal is based on early studies on superstring theory and M-theory, which have sug- 
gested that our four dimensional world is embedded into a higher dimensional spacetime. 
Particularly, the 10 dimensional Eg ® Eg heterotic superstring theory is a low-energy limit 



Bt 



of the 11 dimensional supergravity, under the compactification scheme M^° x S1/Z2 
Thus, the 10 dimensional spacetime is compactified as M^ x CY^ x S1/Z2, implying that 
our Universe (a brane) is embedded into a higher dimensional bulk. In this paradigm, the 
standard model particles are open strings, confined on the braneworld, whilst the gravitons 
and the closed strings can freely propagate into the bulk J5|. 

The Randall- Sundrum Type II model has the virtue of providing a new type of compacti- 
fication of gravity [l|, |2| . Standard 4D gravity can be recovered in the low-energy limit of the 



model, with a 3-brane of positive tension embedded in 5D anti-de Sitter bulk. The covariant 
formulation of the braneworld models has been formulated in 6|, leading to the modification 
of the standard Friedmann equations on the brane. It turns out that the dynamics of the 
early Universe is altered by the quadratic terms in the energy density and by the contri- 
bution of the components of the bulk Weyl tensor, which both give a contribution in the 
energy momentum tensor. This implies a modification of the basic equations describing the 
cosmological and astrophysical dynamics, which has been extensively considered recently 

Q. 



The recent observations of the CMB anisotropy by WMAP 10| have provided convincing 



evidence for the inflationary paradigm ll|, according to which in its very early stages the 



Universe experienced an accelerated (de Sitter) expansionary phase (for recent reviews on 
inflation see [12|). 

At the end of inflation, the Universe is in a cold and low-entropy phase, which is utterly 
different from the present hot high-entropy Universe. Therefore the Universe should be 
reheated, or defrosted, to a high enough temperature, in order to recover the standard Hot 



Big Bang 13|. The reheating process may be envisioned as follows: the energy density in 



zero-momentum mode of the scalar fleld decays into normal particles with decay rate F. The 



decay products then scattered and thermalize to form a plasma 12 1. 

Apart from the behavior of the inflaton field, the evolutions of dark energy and dark 
matter in reheating stage were also considered. In IJ], dark energy and dark matter were 
originated from a scalar field in different stages of the inflation, according to a special form 
of potential. Meanwhile, the conditions for unifying the description of inflation, dark matter 
and dark energy were considered in 15|. A specific model was later proposed in [l6[, by 
using a modified quadratic scalar potential. The candidates of dark matter in [l5| and 16 1 
were oscillations of a scalar field. However, it may be possible that dark matter existed on 
its own without originating from the scalar field. This may pose less stringent constraint on 
the scalar field, so that dark matter can be included in inflation paradigm in a easier way. 
On the other hand, it was proposed that the decay products of scalar field acquired thermal 
mass [17 1 . 

The reheating in the braneworld models has also been considered recently. In the context 
of the braneworld inflation driven by a bulk scalar field, the energy dissipation from the 
bulk scalar field into the matter on the brane was studied in 18|. The obtained results 
supports the idea that the brane inflation model, caused by a bulk scalar field, may be a 
viable alternative scenario of the early Universe. The inflation and reheating in a braneworld 
model derived from Type IIA string theory was studied in [19l|. In this model the inflaton 
can decay into scalar and spinor particles, thus reheating the Universe. A model in which 
high energy brane corrections allow a single scalar field to describe inflation at early epochs 
and quintessence at late times was discussed in 20|]. The reheating mechanism in the model 
originates from Born-Infeld matter, whose energy density mimics cosmological constant at 
very early times and manifests itself as radiation subsequently. The particle production at 
the collision of two domain walls in a 5-dimensional Minkowski spacetime was studied in 



2l| . This may provide the reheating mechanism of an ekpyrotic (or cyclic) brane Universe, 



in which two BPS branes collide and evolve into a hot big bang Universe. The reheating 
temperature Trh in models in which the Universe exits reheating at temperatures in the 



MeV regime was studied in 



221 ■ and a minimum bound on Trh was obtained. The de- 



rived lower bound on the reheating temperature also leads to very stringent bounds on the 
compactification scale in models with n large extra dimensions. The dark matter problem 
in the Randall- Sundrum type II braneworld scenario was discussed in 23|, by assuming 



that the lightest supersymmetric particle is the axino. The axinos can play the role of cold 



dark matter in the Universe, due to the higher reheating temperatures in the braneworld 
model, as compared to the conventional four-dimensional cosmology. The impact of the 
non-conventional brane cosmology on the relic abundance of non-relativistic stable particles 
in high and low reheating scenarios was investigated in [241] . In the case of high reheating 
temperatures, the brane cosmology may enhance the dark matter relic density by many order 
of magnitudes, and a stringent lower bound on the five dimensional scale may be obtained. 
In the non-equilibrium case, the resulting relic density is very small. The curvaton dynamics 



in brane- world cosmologies was studied in 25 1. 



Brane-worlds with non-constant tension, based on the analogy with fluid membranes, 
which exhibit a temperature-dependence according to the empirical law established by 



Eotvos, were introduced in 26|. This new degree of freedom allows for evolving gravi- 
tational and cosmological constants, the latter being a natural candidate for dark energy. 
The covariant dynamics on a brane with variable tension was studied in its full general- 
ity, by considering asymmetrically embedded branes, and allowing for non-standard model 
fields in the 5-dimensional space-time. This formalism was applied for a perfect fluid on 
a Friedmann brane, which is embedded in a 5-dimensional charged Vaidya-Anti de Sitter 
space-time. For cosmological branes a variable brane tension leads to several important 
consequences. A variable brane tension may remove the initial singularity of the Universe, 
since the brane Universe was created at a finite temperature T^ and scale factor a^in J27|. 
Both the brane tension and the 4-dimensional gravitational coupling 'constant' increase with 
the scale factor from zero to asymptotic values. The 4-dimensional cosmological constant is 
dynamical, evolving with a, starting with a huge negative value, passing through zero, and 
finally reaching a small positive value. Such a scale-factor dependent cosmological constant 
has the potential to generate additional attraction at small a (as dark matter does) and late- 
time repulsion at large a (dark energy). The evolution of the brane tension is compensated 
by energy interchange between the brane and the fifth dimension, such that the continu- 
ity equation holds for the cosmological fluid [27|. The resulting cosmology closely mimics 
the standard model at late times, a decelerated phase being followed by an accelerated ex- 
pansion. The energy absorption of the brane drives the 5D space-time towards maximal 
symmetry, thus becoming Anti de Sitter. Other physical and cosmological implications of a 



varying brane tension have been considered in 



It is the purpose of the present paper to further investigate the cosmological implications 



of a varying brane tension. As a first step in our study, we consider a thermodynamic 
interpretation of the varying brane tension models, by showing that the field equations with 
variable A can be interpreted as describing matter creation in a cosmological framework. 
The particle creation rate is determined by the variation rate of the brane tension, as well 
as by the brane-bulk energy-matter transfer rate. In particular, by adopting a theoretical 
model in which the brane tension is a simple function of the scale factor of the Universe, we 
consider the possibility that the early inflationary era in the evolution of the brane Universe 
was driven by a varying brane tension. A varying brane tension may also be responsible for 
the generation of the matter after reheating, as well as for the late time acceleration of the 
Universe. 

The present paper is organized as follows. In Section [ITl we present the field equations of 
the brane world models with varying brane tension and we write down the basic equations 
describing the cosmological dynamics of a flat Friedmann- Robertson- Walker Universe. The 
thermodynamic interpretation of the brane-world models with varying brane tension and 
brane-bulk matter-energy exchange is considered in Section lllli A power-law inflationary 
brane-world model with varying brane tension and non-zero bulk pressure is obtained in 
Section IIVI The analytical behavior of the cosmological model with varying brane tension 
is considered in Section |Vl by using the small and large time approximations for the brane 
tension. The numerical analysis of the model is performed in Section IVIIi We discuss and 
conclude our results in Section IVIIIi 

II. GEOMETRY AND FIELD EQUATIONS IN THE VARIABLE BRANE TEN- 
SION MODELS 

In the present Section we present the field equations for brane world models with varying 
brane tension, and the corresponding cosmological field equations for a flat Robertson- 
Walker space-time. 

A. Gravitational field equations 

We start by considering a five dimensional {5D) spacetime (the bulk), with a large neg- 
ative 5D cosmological constant ^^•'A and a single four-dimensional (4/)) brane, on which 



usual (baryonic) matter and physical fields are confined. The AD braneworld (^^^M, ^^^g^u) 
is located at a hypersurface (-B {X^) = O) in the 5D bulk spacetime (*^^^M, ^^^qab) with 
mirror symmetry, and with coordinates X^, A = 0, 1, ..., 4. The induced AD coordinates on 
the brane are x'^,fi = 0, 1,2,3. We choose normal Gaussian coordinates, and therefore the 
5D metric is related to the AD metric by the relation ^^^Qmn = ^^^ Qmn + ^MnN, where n^ 
is the normal vector. 

The induced AD metric is gu = ^^^gij — n^rij, where n/ is the space-like unit vector field 
normal to the brane hypersurface *^^^M. The basic equations on the brane are obtained by 
projections onto the brane world with Gauss equation, Codazzi equation and Israel junction 
condition, the projected Einstein equation are given by 



where 



^fiu — ^dfiu 'T k 1 nj^ -\- k b^i, €^1, + L + rf^y + r fiu) (J-j 



S,u = \TT,y - \t,^T: + ''-'\^ g,y, (2) 

e^u = CABCDn^n^g^g,, (3) 



and 



F,. =(^) TABg^g^ + i^'^TABU^n^ ~a) ^'"'' 



(4) 



respectively. 

Apart from the terms quadratic in the brane energy-momentum tensor, in the field equa- 
tions on the brane there are two supplementary terms, corresponding to the projection of 
the 5D Weyl tensor e^^, and of the projected tensor F^^,, which contains the bulk matter 
contribution. Both terms induce bulk effects on the brane. 

Also, the possible asymmetric embedding is characterized by the tensor 

L^, = K,,k - K,,K: - ^ (i?2 _ K^pK^^) , (5) 



with trace L = KapK^^ — K^, and trace-free part Z^^f = K^^K — K^^^K^ + Lg^^iy/A, respec- 
tively. 

For a Z2 symmetric embedding Kfj_i, = 0, and thus Z^;^ = 0. P^u is given by the pull-back 
to the brane of the energy-momentum tensor characterizing possible non-standard model 
fields (e. g. scalar, dilaton, moduli, radiation of quantum origin) living in 5D, 

2P7 Z \^^ 

P.. = ^(^>'(^)T./3) , (6) 



which is traceless by definition. Another projection of the 5D sources appears in the brane 
cosmological constant A. which is defined as 

L 2P- 



A = Ao---— (n-n/^(5)T.^), (7) 

where 2Ao = kjX + kjA^. 

In the case of a variable brane tension, the projected gravitational field equations on the 
brane have a similar form to the general case, 

Gfj,p = —Kg^y + k T^y + k S^y — e^y + L^^ + P^y + F^y. (8) 

However, the evolution of the brane tension appears in the Codazzi equation, and in the 
differential Bianchi identity. The Codazzi equation is 



V^K^ - VyK = kl{gPn^i^)Tp^), (9) 

and it gives the conservation equation of the matter on the brane as 

V^T,^ = V.A-A(^X^'^^pa)- (10) 

The differential Bianchi identity, written as V^-Rp^ = |Vp-R, gives 

V^{e,y - L^y - P,y) = ^ + fv.(n^n<^(5)Tp.) - ^^ A^n'T.p) 

+ \{T,yV^T - TVyT)] - |(r,^ - r<)VpA). (11) 

From Eq. Q, one can introduce an effective non-local energy density U, which can be 
obtained by assuming that e^y in the projected Einstein equation behaves as an effective 
radiation fluid, 

- £,,u = YAt/(MpM^ + —V), (12) 

where u^ is the matter four- velocity, and hfj^y = g^j^y + u^y, respectively. 

B. Cosmological models with dynamic brane tension 

We assume that the metric on the brane is given by the flat Robertson- Walker- Friedmann 
metric, with 

^^^g^ydx^dx"" = -dt^ + a\t){dx^ + dy'^ + dz^), (13) 



where a is the scale factor. The matter on the brane is assumed to consist of a perfect 
fluid, with energy density p, and pressure p, respectively. The gravitational field equations, 
governing the evolution of the brane Universe with variable brane tension, in the presence 
of brane-bulk energy transfer, and with a non-zero bulk pressure, are then given by |26|, |27| 
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P^t,^U 



p|l + ^)+3p(l + ^)+2C/ 



0-=' 



fc?^ 



(14) 

(15) 
(16) 
(17) 

(18) 



where P5 describes the bulk-brane matter-energy transfer, while Pb is the bulk pressure. 

An important observational parameter, which is an indicator of the rate of expansion of 
the Universe, is the deceleration parameter q, defined as 



Q 



d 
di 



I 



aa 



a? 



d/a 

[d/ay 



(19) 



If g < 0, the expansion of the Universe is accelerating, while g > indicates a decelerating 
phase. 

III. THERMODYNAMIC INTERPRETATION OF THE VARYING TENSION IN 
BRANE-WORLD MODELS 



For the sake of generality we also assume that there is an effective energy-matter trans- 
fer between the brane and the bulk, and the brane-bulk matter-energy exchange can be 
described as 

abb [(^OV"" rr ron\ 

P5 = — —per (^— j H, (20) 

where abb is a constant, pcr is the present day critical density of the Universe, and oq is the 
present day value of the scale factor. 

In the presence of a varying brane tension and of the bulk-brane matter and energy 



exchange, the energy conservation equation on the brane can be written as 



p + ?,{p + p)H 



-p\ UbbH 



(21) 



where we have used Eq. (1201) for the description of brane bulk energy transfer, by taking into 
account that p = pcr (ao/a) "". We suppose that the matter content of the early Universe 
is formed from m non-interacting comoving relativistic fluids with energy densities and 
thermodynamic pressures pi(t) and Pi{t), i = 1, 2, ..., m, respectively, with each fluid formed 
from particles having a particle number density nj(t), i = 1,2, ..., m, and obeying equations 
of state of the form pi{t) = kiuj' , pi{t) = {ji — I) Pi, i = 1,2, ...,m, where ki = pQi/nl\ > 
0, i = 1,2,.. .,m, are constants, and 1 < '~fi < 2, i = 1,2,. ..,m. For example, we can 
consider that the particle content of the early Universe is determined by pure radiation (i. 
e., different types of massless particles, or massive matter (baryonic and dark) in equilibrium 
with electromagnetic radiation and decoupled massive particles. The total energy density 
and pressure of the cosmological fluid results from summing the contribution of the / simple 
fluid components, and are given by p{t) = Y^=iPi{'^) ^^^1 p{t) = Yl^iPii't): respectively. 
For a multicomponent comoving cosmological fluid and in the presence of variable brane 
tension and bulk-brane energy exchange, Eq. fl2T]) becomes 



5^[A + 3(a+p,)^] = -5Zp.W 



A 



ET=lP^it) 



- CibbH 



Eq. 



can be recast into the form of m particle balance equations, 
hi{t) + 3ni{t)H = Ti{t)ni{t),i = 1, 2, ..., m. 



(22) 



(23) 



where Ti{t), i = 1,2, ..., m, are the particle production rates, given by 



r.(t) = -- 

7i 



A 



OihhH 



1,2, ...,m. 



(24) 



mpiit) 

In order for Eq. (123|) to describe particle production the condition Vi{t) >Q,i = 1,2, ...,m, 
is required to be satisfied, leading to the following restriction imposed to the time variation 
rate of the brane tension 

X<abblpi{t)H,i = 1,2,.. .,m. (25) 

Note that if Ti{t) = Q, i = 1,2, ...,m, we obtain the usual particle conservation law of the 
standard cosmology. Of course, the casting of Eq. (I22l) is not unique. In Eqs. (I23l) and (l2ll) . 



9 



we consider the simultaneous creation of a multicomponent comoving cosmological fluid, 
but other possibilities can be formulated in the same way (for example, creation of a single 
component in a mixture of fluids). 

The entropy Si, generated during particle creation at temperatures Tj, i = 1, 2, ..., -m, can 
be obtained from Eq. f l23p . and for each species of particles has the expression 



T,;. 



dSi 
If 



1 

7i 



A 



mpi{t) 



- OihbH 



Pi{t)V,i = 1,2, ...,m. 



(26) 



where V is the volume of the Universe, or, equivalently, 

dSi liPi{t)V 



-Ti{t),i = l,2,...,m. 



(27) 



dt Ti 

In a cosmological fluid where the density and pressure are functions of the temperature 
only, p = p (T), p = p (T), the entropy of the fluid is given by 5 = {p + p) V/T = •yp (t) V/T. 
Therefore we can express the total entropy S{t) of the multicomponent cosmological fluid 
filled brane Universe as a function of the particle production rate only, 

ft 



Sit) = J2So^exp / T,it')dt' 



(28) 



where Soi > 0, i = 1,2, ...,m, are constants of integration. In the case of a general perfect 
comoving multicomponent cosmological fluid with two essential thermodynamical variables, 
the particle number densities rii, i = 1, 2, ..., m, and the temperatures Tj, i = 1,2, ..., m, it is 
conventional to express pi and pi in terms of Ui and Tj by means of the equilibrium equations 
of state Pi = Pi {rii, Ti), pi = pi {rii, Ti), i = 1, 2, ..., m. By using the general thermodynamic 

relation 

dpi Pi + Pi Ti dpi 



A = 1,2, ...,m. 



drii Ui Hi dTi ' 

in the case of a general comoving multicomponent cosmological fluid Eq. f l2T|) can also be 
rewritten in the form of m particle balance equations. 



hiit) + 3ni{t)H = ri(t)ni{t),i = 1, 2, ..., m. 



(30) 



with the particle production rates Ti{t) given by some complicated functions of the thermo- 
dynamical parameters, brane tension and brane-bulk energy exchange rate. 



r.(t) 



pi 



pi+Pi 



A d In Pi I Ti 2 ^i 

- — -r - abbH + Ti \ — - Ci — 

lpi{t) dTi \Ti Hi 



,i = 1,2, ...,m. 



(31) 
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where Cf = (dpi/dTi) / (dpi/dTi). The requirement that the particle balance equation 
Eq. fl30|) describes particle production, Ti(t) > 0, i = 1,2, ...,m, imposes in this case the 
following constraint on the time variation of the brane tension, 

A < a,M{t)H-lMm^ i^-C^^] ,^ = l,2,...,m. (32) 

oTi \Ti Hi J 

In the general case the entropy generated during the reheating period due to the variation 
of the brane tension and the bulk-brane energy exchange can be obtained for each component 
of the cosmological fluid from the equations 



dSj _ {pi+Pi)V 
dt ~ Ti 



A din Pi [ Ti 2 'h' 



abbH + Ti—-f^ 777 -G 



2 "-« 



1,2,. ..,m, (33) 



lpi{t) OTi \Ti ' rii 

while the total entropy of the Universe is given by S{t) = J2iLi Siit). 

The entropy flux vector of the kth component of the cosmological fluid is given by 

5W" = nfcafcM",A; = l,2,..,m, (34) 

where ak, k = 1,2, ...,m, is the specific entropy (per particle) of the corresponding cosmo- 
logical fluid component and u" is the four- velocity of the fluid. By using the Gibbs equation 
nTda = dp— [{p + p) /n] dn for each component of the fluid, and assuming that the entropy 
density a does not depend on the brane tension, we obtain 

S\t^- = -;i (a - abbHp) - ^^^, k = l,2, ..,m, (35) 

where /i^ is the chemical potential deflned by /i^ = [{pk + Pk) /nk\ ~ Ti^ai^. The chemical 
potential is zero for radiation. For each component of the cosmological fluid the second law 
of thermodynamics requires that the condition 

^g)">0,fc = l,2,..,m, (36) 

has to be satisfled. 

IV. POWER LAW INFLATION IN BRANE WORLD MODELS WITH VARYING 
BRANE TENSION AND BULK PRESSURE 

For a vacuum Universe with p = p = 0, in the presence of a non-zero bulk pressure and 
matter- energy exchange between the brane and the bulk, the fleld equations Eqs. (JH]) take 

11 



the form 



and 



1 (Pa P' . , \ 

^-T^ = t; -PB-lu, (38) 

a dr^ 2 



'' (lua')^-a'-^(^!^-pBy (40) 



(ir 
where we have introduced a set of dimensionless variables {t,1,pb,u,P5) defined as 

/3 3 3 3 3 

Moreover, we consider that the five- dimensional cosmological constant A5 = 0. We assume 
that the inflationary evolution is of the power law type, and therefore a = r", where a is a 
constant. Then Eqs. (!37|) and (138!) give 



'/2 \ 3a (2a - 1 

and 



n--VB\=-^^^^^iP^. (42) 



llu = ^. (43) 



r2 



Eq. ( Hoi) is identically satisfied. In order to completely solve the problem, we need to 
specify the form of the energy matter-transfer from the bulk to the brane. By assuming a 
functional form given by ps = PqqT^^, where /3 > and po5 > are constants, we obtain 
immediately 

/ (r) = ,llj—r-^^+^ „^ (r) = 1^ r^^C/^-i) - ^^ ^^^ " ^) 

^' V 3/3-1 '^""^^ 3(/3-l)2^ 2r2 ' 

/27 
n{r) = ^-a{P-l)r^-'. (44) 

The Hubble parameter of the Universe during the inflationary phase is given by iJ = a/t. 
The deceleration parameter is obtained as g = rf {1/H) /dt — 1 = (1 — a) /a. Therefore, if 
a > 1, g < 0, and the brane world Universe experiences an inflationary expansion. 



12 



V. SCALE FACTOR DEPENDENT BRANE TENSION MODELS 

In the following we assume that there is no matter-energy exchange between the bulk 
and the brane, P5 = 0, and that the bulk pressure is also zero, Pb = 0. For the matter on 
the brane we adopt as equation of state a linear barotropic relation between density and 
pressure, given by 

p={w-l)p, (45) 

where w = constant and w G [1,2]. Therefore Eq. flT6|l gives 

p + 3Hwp = -A, (46) 

while Eq. flTTl) gives immediately 

XU = ^, (47) 

where f/o is an arbitrary constant of integration. In the following, in order to simplify the 
analysis, we assume that Uq = 0. 

In order to explain the main observational features of modern cosmology (inflation, re- 
heating, deceleration period and late time acceleration, respectively), we assume that the 
brane tension varies as a function of the scale factor a according to the equation 

A^ = Xle-'^^' - V^ + A?, (48) 

where /3, Aq and Ai are constants. 

Suppose tin, flin and pin are the values of the time, of the scale factor, and of the energy 
density before inflation. Generally, in the present paper we use the subscript "in" to denote 
the values of the cosmological parameters before the inflation, and the subscript "en" to 
denote values after inflation. Thus, for example, iV = In (aen/ffin) is the e-folding number. 
The basic physical parameters of our model are tin, Qin, Pin, N,k5, ^^^A, Aq, Ai, and (3, 
respectively. The coupling constant k^ and the five-dimensional cosmological constant '^^^A 
are constrained by the present value of the gravitational constant. 



6 V kl ' "' "^V 6 



kl /_ 6(5)A ^ ^2 ^ kl\ -^ = SttG ^ 1.68 X 10-55eV-2, (49) 



and by the constraints on the 5D cosmological constant J9| 



I9 



^(5)A ^ ^—— ^ -2.3 X lO-^eV^, (50) 

2 (0.1mm)2 ' ^ ^ 
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where we have used the natural system of units with h= c = 1. /^From these two conditions, 
we obtain fc| ^ 3.6 x 10"^°^ eV"^ and (^^A ^ -7.7 x 10^^ eV^ Besides, the value of Ai can 
be obtained from the value of the present day dark energy pdark ~ lO^^^eV^ (8|, 

^A? = STTGpdark ~ 1.6 X lO^^W^, (51) 

which gives Ai ~ 1.6 x 10^^ eV^. We can have a backward checking on Eq. (HOl) . from which 
it follows that the condition A^ ^ — 6'^^^A/A;| is indeed satisfied. We also choose Aq to be of 



the same order of magnitude as the vacuum energy p^ac "^ 10^°° eV^ at GUT scale [8|, |33| . 



Ma2 = STrGp^ae ^ lO^^eV^, (52) 



which gives Aq ^ 1.3 x 10^^ eV. The scale difference between Aq and a/— 6(^)A/A;| + Af is 



Aq/a/— 6(^)A//i;| ^ 10^^. The differences in the scales of Aq and Ai are of the order of ^ 56. 
For the remaining model parameters tjn, flm, Pin? N, /3, we constraint them in the next section. 
When a is very small, the brane tension A ~ Aq dominates the early Universe at the time 
of inflation. Due to the exponential expansion of the Universe, the brane tension quickly 
decays to a constant just after the inflation. The decay of the brane tension will generate 
the matter content of the Universe, according to Eq. ( TT6l) . This happens also during the 
accelerated expansion period of the Universe. Matter is created during all periods of the 
expansion of the Universe, but the most important epoch for matter creation is near the 
end of inflation. In the evolution of the Universe there is one moment when a = 0, which 
corresponds to the moment when the Universe switches from the accelerating expansion to 
a decelerating phase. After the matter (which is mainly in the form of radiation) energy 
density reaches its maximum, the Universe enters into a radiation dominated phase, and the 
quadratic term in Eq. flT^ will become dominant first. The matter energy density continue 
to decrease due to the expansion. When the linear term in matter equals the quadratic 
term, the Universe switches back to the ACDM model. Therefore, the Universe enters in 



the matter dominated epoch at about 4.7 x lO^yr 3J]. Then the matter term equals the 



residue term in Eq. (8) at about lOGyr [34]. This is the second moment in the evolution 
of the Universe when a = 0. After this moment, the Universe enters in an accelerating 
expansionary phase again, and its dynamics is controlled by the term Ai. 
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VI. QUALITATIVE ANALYSIS OF THE MODEL 

In the present Section we consider the approximate behavior of the cosmological model 
with varying brane tension in the different cosmological epochs. 

A. Early inflationary phase: 2/3a^ <^ 1 

When the scale factor a is very small, the exponential factor e"^^*^ in Eq. (1481) can 
be approximated by 1. Therefore, the brane tension is given by A^ ~ Aq, and physically it 
corresponds to the vacuum energy necessary to give an exponential inflation. Since /c^Aq/G ^ 
^^^A, from Eq. (TT^ the scale factor evolves in time as an exponential function of time, given 
by 

a = aine('^i^°/6)*, (53) 

where a^^ is the value of the scale factor prior to inflation. The e-folding number is given 
by A^ = ln(a/ain), which should be roughly of the order of A^ > 60 in order to solve the 
flatness. Horizon problem, etc. In the present paper we adopt for A^ the value A^ = 70. Since 
2(3al^ ^ 1 at the end of inflation, we can roughly estimate the value of Gen to be 

1 



flpTl ^ 



From the value of A^ we adopted, we obtain the value of a^^ as 



(54) 



According to Eq. f l53|) . the end time of the inflation ten can be estimated as 

A;5Ao(tcn - tin)/6 ^ klXot.jQ ^ N, (56) 

which implies that ten ~ 10~^^ s. The value of tin is insensitive to the variation of the initial 
conditions, provided that tin is at least one order smaller than ten- With the adopted value 
of the e-folding A^ and /3, we can fix the values of a^^ and of tin, respectively. Since at the 
beginning of the inflationary stage the matter is not yet generated, we have pin = 0. For the 
deceleration parameter, from Eq. ( fT9l) we find 
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-1. (57) 



By substituting the tension we can rewrite Eq. (H6|) as 

dp ^ Ida dX 2/3A^aae-2/3a2 

-— + 3w — rP= — r = — , ^- (58) 

dt a dt dt /\2„-2/3a2 6(5) A I \2 

For 2/3a^ ^ 1, the exponential factor does not change much as a increase. Therefore, 

^ + 3^i^p = , '^^°"' - ^^ = ^klXy/?,. (59) 



^=5 

where we have also used Eq. fl53|) . Conversion of the matter from the brane tension energy 
begins already during the inflationary stage, and the rate of the conversion is proportional to 
a? during inflation. Therefore, the matter density also rises exponentially in the late stages 
of the inflationary phase. The matter generation rate becomes most important at the end 
of inflation. 

B. Reheating period: 2/3a^ !=s 1 

During the reheating period the evolution of the matter gradually changes from an ex- 
ponential increase to a power law decrease, p oc a^^'" . In our model the matter density 
is a smooth function, which is strictly increasing from the beginning of inflation, and then 
strictly decreases after the end of the reheating phase. Therefore there must be a maximum 
value of the density pmax at a time tmax- After tmax, the Universe was dominated by matter 
which is in the form of radiation, and almost all the energy of the brane tension converted 
into matter. The temperature of matter, corresponding to a radiation dominated Universe 
at tmax, is denoted Trh, and is given by 

Pmax = — (^^Rh)^, (60) 



where k is Boltzmann's constant. Current theory on gravitinos production constraints Trh 



to be Tt,h < 10^ - 10^° GeV 



30 



32| . Therefore the maximum density of the Universe must 
satisfy the condition pmax < 10'^^ — 10^° GeV^. The maximum density can be obtained from 
the condition p\t=t^^^ = 0, and, with the use of Eq. (14) it is given as a solution of the 
equation 



Ida _ dX _ 2(3Xlade-^^"'^ 

O iV ' /'max — 7 — , '■ 

a dt dt /\2„-2/3a2 6(5) A , \2 



(61) 



t=tn 
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With the rough approximation /3a^|f=f^^,^ ~ (3al^ ^ 1, Eq. ( 161]) can be written as 



Swpn 



^Age-2/3a^ - ^ + A? 



2A2e-2/^«' 



t=t„ 



^2g-2/3a^ - ^ + A? 



0.74Ao. (62) 



t=t„ 



This relation gives the maximum matter density of the Universe. And the value of Aq ^ 
lO^^GeV^ that we have chosen is consistent with the maximum value of the matter energy 
density. 

C. Matter Domination period: 2/3a^ >> 1 and 2Ap » Xf 

At this stage, the Universe is dominated by matter and the brane tension is roughly a 
constant. During this period, the key difference with the conventional cosmological models 
is the presence of the quadratic term in matter density, which will dominate the dynamics 
of the Universe at the beginning of this period. The evolution equation of the scale factor is 

I . 4 (V2A?T7) « 4p, (63) 

and the evolution of the matter density is given by 

dp Ida 
dt a dt 

where we have assumed that immediately after the matter energy reaches its maximum the 
matter is in the form of radiation. The solution is p = constant/a'^ = PmaxO.'^(imax)/a^ ~ 
Ao//3^a^. Therefore, the scale factor and deceleration parameter evolve as 

ait) = ( ^Pmaxa'(tmax)t ) ^ ~ ^t, (65) 



I 6p 



ad %f rp(l + ^) + 3(w;-l)p(l + f;, , , 

q = -— = ^^^^ 74^. — ~ 1 + 3w = 5. (66) 



IS VA^ ^ A 



18 ^^ ' \ ' 

2 



After the quadratic density and linear density equality 2Ap = p , the linear term in matter 
will take over, i.e. 2Ap ^ p^. The Universe enters in the ACDM model at this radiation 
dominated phase, and its dynamics is described by the equations 

(67) 
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and 



dp Ida , . 

dt a dt 



respectively. The solution for the scale factor is 

with p = PmaxO^(i^max)/'^'^- During this period the deceleration parameter evolves as 



fc|A 
36 



[p + 3{w - l)p] 3w-2 



After this stage, the Universe switches from radiation dominated to non-relativistic matter 

n 

dominated at trm = 1.5 x 10^^ s [34]. To simulate the transition from the radiation dominated 
period to the baryonic matter dominated period one can introduce a time varying w given 
by y 

4trm/3 + t 

w = . (71) 

f -\-f 

In the matter dominated era {w = 1), the deceleration parameter shifts to q = 0.5, 
according to Eq. ( 1701) . Recall that /3 is still free, and we can use it to match the scale factor 
at the radiation-matter equality. With the use of Eq. ( l69l) . and by taking arm = 2.8 x 10~^ 
J34| . we obtain 



1/2 
2.6 X lO-'^'s-' ^ I ■^'\:,r'' 1 , (72) 



18/32 



which gives Aq ~ 10 -"^^ x /3^. This gives the estimate of /3 ^ 10^^. 

D. Dark Energy Domination era (2Ap < Af ) 

With the increase of the cosmological time, the constant term Ai in Eq. (HM will dominate 
over matter. Thus this term plays the role of the dark energy of the standard ACDM models. 
In its late stages of evolution the Universe becomes "dark energy" dominated. The Universe 
turn to an exponential acceleration again, with 

aoce(^i/6)^i*, (73) 

but with a time scale much longer than the inflationary time scale. The deceleration pa- 
rameter will converge to 



showing that the expansion of the universe is accelerating. 

VII. NUMERICAL ANALYSIS OF THE MODEL 

The field equation can be rewritten in a simple form by introducing as set of dimensionless 
variables r, /, and r, defined as 



r = h^^-^t,p = A;5-V3x(-A5)r,A = A;^- VS x (-A5)/, (75) 

respectively. Then the field equations Eqs. (TT^ - (TTS]) can be written in a dimensionless form 

as 



Id^a 1 P 1. r /. 2r\ _ .. /. r 



a dr'^ 3 6 6 



r\l + —j+3{w-l)r(^l , 



(77) 



dr 1 da dl , , 

— + 3w-—r = —-. 78) 

dr a dr dr 

By rescaling the four- dimensional cosmological constant so that ^^•'A = [kl{—A^)/2]li.ff, it 

follows that les = — 1 + /^/2. In the dimensionless variables, the deceleration parameter is 

given by g = — {ad'^a/dr^) / {da/dr) , and can be explicitly expressed as a function of the 

physical parameters of the model as 

_ 1 - /V2 + / [r (1 + 2r/l) + 3{w - l)r (1 + r/l)] /2 
^~ -l + P/2 + lr + ry2 ' ^ ^ 

To obtain the numerical solution, one should also consider the rescaled form of A, 

l^ = 2{lle-^''''' + I + II), (80) 

where /q and h are constants, with Zq ~ 10^^ ^ /^ ^ lO^'^^. Eq. ( 17T]) becomes 

4r,m/3 + r 

with Trm = (4.8 X 10^'^ eV)trm = 10^^. Notc that the conversion in Eq. ( 1751) can be written 
out numerically 

t = (208 eV"^)r = (1.4 x 10"^^ s)r, p = (2.0 x 10^° eV)r, A = (2.0 x 10^° eV)l. (82) 

The time variations of the scale factor, of the energy density, of the brane tension, and 
of the deceleration parameter of the Universe are presented, for different scales of vacuum 
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energy Iq, in Figs. [T] respectively. From the plot of a and q, we find that there are 5 
stages in the evolution of the Universe. Namely, the inflation and reheating stage, the 
quadratic density stage, the radiation domination stage, the non-relativistic matter stage, 
and the late time acceleration stage. Comparing plots with different Zq's, we see the effects 
of different vacuum energy scales on the evolution of the Universe. According to the plot 
with different /q'Sj "we find that a larger vacuum energy would provide a faster inflation, and 
it also generates more matter. The matter energy density reaches its maximum at earlier 
times. Although there are many changes in the evolution of the Universe caused by changing 
Iq, these characteristics are hard to be constraint by observations. On the other hand, /q 
cannot be arbitrary due to the theoretical constraint, e.g. vacuum energy density, gravitinos 
production, etc. 





FIG. 1: Time variations of the scale factor a, deceleration parameter q, and energy density r of 
the Universe in braneworld models with varying scale factor dependent brane tension. Different /q 
are plotted: Iq = 10^^ (dotted curve), Iq = 10^^ (solid curve), and Zq = 10^^ (dashed curve). 

The time variations of the scale factor, and of the deceleration parameter of the Universe 
are presented, for different values oi (3, in Figs. [2] respectively. If we assume Aq is fixed 
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by the vacuum energy scale, the value of /3 is well confined. Varying /3 would afi^ect the 
observational constraint of a. From the graph of r, we can cross check that the value of 
Pmax indeed fulfills the requirement of Eq. ( p7]) . Besides /3, we could examine the effect of 
adopting a different e-folding A^. Since a^n is determined by /3 according to the condition 
Eq. f l54p . it follows that it is determined by observational constraints. Different e- foldings 
could be a result of differences in Oin. However, we find from Fig. [3] that a\n does not affect 
the post-infiationary epochs. Therefore Oin is not a robust parameters, and we cannot fully 
determine it. 





1.0x10"'^ 1.2x10"^' 1.4x10" 



1.8x10 2.0x10 



FIG. 2: Time variations of the scale factor a, deceleration parameter q, and energy density r of 
the Universe in braneworld models with varying scale factor dependent brane tension. Different (3 
are plotted: j3 = 10^'^ (dotted curve), j3 = 10^^ (solid curve), and /3 = 10^'^ (dashed curve). 



VIII. DISCUSSIONS AND FINAL REMARKS 

In the present paper we have considered some cosmological implications of the braneworld 
models with variable tension. We have considered a thermodynamic interpretation of the 
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FIG. 3: Time variations of the scale factor a of the Universe in braneworld models with varying scale 
factor dependent brane tension. Different Om are plotted: Om = 10^^^ (dotted curve), Cm = 10^^^ 
(solid curve), and a^i = 10^^^ (dashed curve). 

model, and we have shown that from a thermodynamic point of view a variable brane 
tension can describe particle production processes in the early Universe, as well as the 
entropy production in the early stages of the cosmological evolution. A simple power-law 
inflationary cosmological model has also been obtained. By adopting a simple analytical 
form for A, we have obtained a complete description of the dynamics and evolution of the 
Universe from the stage of the inflation to the phase of the late acceleration. Moreover, the 
differences between the energy scales of the theoretical vacuum energy at inflationary epoch 
and during late time acceleration are studied. 

A variable brane tension can drive the inflationary evolution of the Universe, and also be 
responsible for matter creation in the post-inflationary phase (the reheating of the Universe). 
In the model we have adopted the brane tension converges to a constant that gives the 
gravitational constant, and the residue after cancelation with the A5 term would give the 
dark energy. It is interesting to note that, as opposed to the standard cosmological scenarios, 
in the present model matter creation takes place during the entire inflationary phase, but 
reaches a maximum after the exponential expansion of the Universe ends. Therefore in the 
braneworld models with varying brane tension there can be no clear distinction between 
inflation and reheating. By using numerical analysis as well as approximate analytical 
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methods, we have obtained the result that in this variable tension model there are 5 phases 
in the cosmological evolution of the Universe. At the beginning, there are the inflationary 
and the reheating phases. During these phases, the brane tension is the dominant energy 
in the Universe, and its magnitude is of the order of the vacuum energy at GUT scale. As 
the brane tension decays, the matter is created, and the Universe enters in the hot Universe 
phase of the Big Bang picture. The third phase is the quadratic density domination phase. 
This is a unique characteristic of the brane world scenario, during which the Universe is 
dominated by a quadratic term in the energy density of the radiation. The deceleration 
parameter increases from —1 during inflation to 5 at this phase. To study the cosmological 
dynamics we have introduced a set of dimensionless quantities, which describe the evolution 
of the scaled densities in terms of a dimensionless time parameter r. By using the results 
of the numerical simulations for the rescaled variables, one can obtain some constraints on 
the physical parameters of the model. 
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